QCD prediction for heavy boson transverse momentum distributions 



Jianwei Qiu and Xiaofei Zhang 
Department of Physics and Astronomy, Iowa State University; Ames, Iowa 50011, USA 

(December 4, 2000) 



O 
O 
O 

(N 
O 

Q 

in 



> : 

00 . 
in ■ 
o ■ 
<n : 

o ■ 
o ■ 

i 

CD . 
43 ■ 



X 



We investigate the predictive power of Collins, Soper, and 
Sterman's 6-space QCD resummation formalism for trans- 
verse momentum (Qt) distributions of heavy boson produc- 
tion in hadronic collisions. We show that the predictive power 
of the resummation formalism has a strong dependence on the 
collision energy %/S in addition to its well-known Q 2 depen- 
dence, and the vS dependence improves the predictive power 
at collider energies. We demonstrate that at Tevatron and the 
LHC energies, the Qt distributions derived from 6-space re- 
summation are not sensitive to the nonperturbative input at 
large b, and give good descriptions of the Qt distributions of 
heavy boson production at all transverse momenta Qt < Q- 
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With new data from Fermilab Run II on the horizon 
and the LHC in the near future, we expect to test Quan- 
tum Chromodynamics (QCD) to a new level of accuracy, 
and also expect that a better understanding of QCD will 
underpin precision tests of the Electroweak interactions 
and particle searches beyond the Standard Model |jj . As 
pointed out in Ref. the description of vector and 
scalar boson production properties, in particular their 
transverse momentum (Qt) distribution, is likely to be 
one of the most intensively investigated topics at both 
Fermilab and the LHC, especially in the context of Higgs 
searches. The main purpose of this Letter is to quantita- 
tively demonstrate the predictive power of QCD resum- 
mation formalism for the Qt distributions of heavy boson 
production at Tevatron and the LHC. In particular, we 
concentrate on the small transverse momentum region: 
Qt < Q, where the bulk of the data is. This region is 
also most relevant to the hadronic Higgs production. 

When Qt <C Q, the Qt distribution calculated in con- 
ventional fixed-order perturbation theory receives large 
logarithm, ln(Q 2 / Q 2 -,) , at every power of a s , which is a 
direct consequence of emissions of soft and collinear glu- 
ons by incoming partons. Therefore, at sufficiently small 
Qt, the convergence of conventional perturbative expan- 
sion in powers of a s is impaired, and the logarithms must 
be resummed. 

Resummation of the large logarithms can be carried 
out either in Q^-space directly, or in the impact parame- 
ter, 6-space, which is the Fourier conjugate of Qr-space. 
It was first shown by Dokshitzer, Diakonov, and Troyan 
that in the double leading logarithm approximation, the 
dominant contributions in small Qt region can be re- 
summed into a Sudakov form factor g]. By imposing 
the transverse momentum conservation without assum- 
ing the strong ordering in transverse momenta of radiat- 



ing gluons, Parisi and Petronzio introduced the 6-space 
resummation method which allows to resum some sub- 
leading logarithms J3|. Using the renormalization group 
equation technique, Collins and Soper improved the 6- 
space resummation to resum all logarithms as singular 
as, \^ m (Q 2 /Q 2 T )/Q 



as Qt — * 14] . Using this renor- 
malization group improved 6-space resummation, Collins, 
Soper, and Sterman (CSS) derived a formalism for the 
transverse momentum distributions of vector boson pro- 
duction in hadronic collisions ||. This formalism, which 
is often called CSS formalism, can be also applied to the 
hadronic production of Higgs bosons jjj . 

For Drell-Yan vector boson (V — r )*,W ± ,Z) produc- 
tion in hadronic collisions between hadrons A and B, the 
CSS formalism has the following generic form H 



dcFA+B- 



■v+x 



dQ 2 dydQ 2 T 



(2tt) 2 
Y(Qt, 



d 2 be lQTb W(b,Q,x A ,x B ) 
,x A ,x B ), (1) 



where W term gives the dominant contribution when 
Qt "C Q, while Y term gives corrections that are negligi- 
ble for small Qt, but become important when Qt ~ Q- 
In Eq. (0), x A = e y Q/VS and x B = e - * Q/VS with the 
rapidity y and collision energy y/S. The W in Eq. (|l|) 
includes all powers of large logarithms from ln(l/6 2 ) to 
ln(Q 2 ) and has the following form || 



W(b, 



,Xa,Xb) 



-S(b,Q)w(b,c/b,x A ,x B ) 



(2) 



where c is : 
5(6, Q) = I 



one 



constant of order 



>/b* 7? 



|,§, and 

B(a s (fj,)) , 

with A(a s ) and B(a s ) perturbatively calculable ||. The 
W(b, c/b, xa, xb) in Eq. (Q) depends on only one momen- 
tum scale, 1/6, and is perturbatively calculable as long as 
1/6 is large enough. The large logarithms from ln(c 2 /6 2 ) 
to ln(Q 2 ) in W(b, Q, xa, xb) are completely resummed 
into the exponential factor exp[— S(b, Q)]. 

Since the perturbatively resummed W(b, Q, xa, x b ) in 
Eq. (^|) is only reliable for the small 6 region, an ex- 
trapolation to the large 6 region is necessary in order to 
complete the Fourier transform in Eq. ([!]). In the CSS 
formalism, a variable 6* and a nonperturbative function 
F NP (b, Q ,xa,xb) were introduced ||, 



W C ss{b,Q,x A ,x B ) 



= W(K,Q,x A ,x B ) 
x F NP {b,Q,x A ,x B ) 



where 6* = b/ + (b/b max ) 2 < b max = 
and F NP has a Gaussian-like form, F NP 



(3) 

0.5 GeV" 1 , 
^ exp(— nb 2 ) 
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and the parameter k has some dependence on Q 2 , xa, 
and xb- The data are not inconsistent with such a form 
Jt| — |l0|] - However, improvements are definitely needed for 
the precision tests of the theory |p],pT[ . 

Although the 6-space resummation formalism has been 
successful in interpreting existing data, it was argued 
Jj|,[il"l that the formalism has many drawbacks associ- 
ated with working in impact parameter space. As listed 
in Ref. Q , the first is the difficulty of matching the re- 
summed and fixed-order predictions; and the second is to 
know quantitative difference between the prediction and 
the fitting because of the introduction of a nonperturba- 
tive F NP . In viewing of these difficulties, major efforts 
have been devoted to resum the large logarithms directly 
in Q T -space @Jll[. 

In the following, we argue and demonstrate that both 
of these drawbacks can be overcomed. We show that 
the 6-space formalism works smoothly for all Qt < Q- 
We demonstrate that the Q 2 and V~S dependence of the 
resummed 6-space distribution ensure that the Fourier 
transform is completely dominated by the small 6 region 
in high energy collisions, and consequently, the Qt dis- 
tribution is insensitive to the details of F NP . 

It was known jl| that the 6-space resummed Qt dis- 
tribution from Eq. (|l|) becomes unphysical or negative 
when Qt is large. For example, a matching between the 
resummed and fixed-order calculations has to take place 
at Qt ~ 50 GeV for W production when these two pre- 
dictions cross over pJ| . On the other hand, we expect 
that the predictions given by the 6-space resummation 
in Eq. (|l]) should work better when Qt is large, because 
the perturbatively calculated Y term dominates and the 
predictions should be less sensitive to W and its non- 
perturbative input. We find that this puzzle was mainly 
caused by the lack of numerical accuracy of the Fourier 
transform from the 6-space to Qr-space. 

Since there is no preferred transverse direction, the 
W in Eq. (@) 

is a function of 6 = |6|, and the Fourier 
transform can be written as 



d 2 be iQT b W{b,Q,x A ,x B ) 



(2n) 2 



1 f°° r 
- dbbJo(Q T b)e~ s< - b ^W(b,-,x A ,x B ), (4) 



2tt 



where Jo(z) with z = Qt6 is the Bessel function. Because 
of the oscillatory nature of the Bessel function, high ac- 
curacy of the numerical integration over 6 is crucial for 
a reliable result. The number of oscillations strongly de- 
pends on the value of Qt for the same range of 6. For 
example, when 6 £ (0, 2) GeV -1 , JoiQrb) crosses zero 0, 
6, and 63 times for Qt — 1, 10, and 100 GeV, respec- 
tively. It is clear that numerical accuracy is extremely 
important for the large Qt region. We noticed that 
most work published in the literature used some kind 
of asymptotic form to approximate the Bessel function 
when z = Qrb is large. We believe that the asymp- 
totic form is a source of the uncertainties observed for 



the large Qt region. Instead of using an asymptotic 
form, we use an integral form for the Bessel function 
Jq(z) — ~ Jq cos {z sin(0)) d6. The great advantage of 
using an integral form is that we can control the nu- 
merical accuracy of the Bessel function by improving the 
accuracy of the integration. With the integral form of 
the Bessel function, we show below that the 6-space re- 
summed Qt distributions are smoothly consistent with 
data for all transverse momenta up to Q ||. 

The predictive power of the CSS resummation formal- 
ism relies on the fact that the integration over 6 in Eq. (Q) 
is dominated by the region where 6 ~ 1/Q, because the 
exponential factor exp[— 5(6, Q)} in Eq. (||) suppresses the 
6-integral when 6 is larger than 1/Q Using the saddle 
point method, it was shown |||| that even at Qt = 0, 
the 6-integration in Eq. (^) is dominated by an impact 
parameter of order 



9SP 



A 



QCD 



A 
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where A = 16/(49 — 2n/) « 0.41 for quark flavors n/ = 5. 
When Q ~ Mz or Myy, the momentum scale c/bsp in 
Eq. (|4|) is of a few GeV, at which the perturbation theory 
is expected to work. Consequently, the predictive power 
of the CSS formalism should not be very sensitive to the 
F NP at large 6 as long as Q is large. 

The 6sp in Eq. (^J) was derived with an assumption that 
the 6-dependence in W(b, c/b, xa, xb) is smooth around 
6sp. We find that the 6-dependence in W(b, c/b, xa, xb) 
is important. When xa and xb are small, this 6- 
dependence reduces the numerical value of saddle point, 
and consequently, increases the predictive power of the 
6-space resummation formalism 0]. 

Taking into account the full 6-dependence of 
W(b,c/b,XA,XB), the saddle point for the 6-integration 
in Eq. (|J) at Qt — is determined by solving the follow- 
ing equation 



In (6e~ s ^) + - In (w(b, ^,x A ,x B )) = 0. (6) 



d 

Jb*"V~ ) ' db"' V" v "' 6 

The 6sp in Eq. (||) corresponds to a solution of Eq. (^) 
without the second term and keeps only the first order 
term of the A(a s ) in S(b,Q). As shown in Ref. ||, the 
second term in Eq. (|J) is proportional to the evolution of 
parton distribution <fi(x,[i): — jj^jt) < M a; ' i l = §)• The 
evolution (d/d\n fi)(j)(x, fi) is positive (or negative) for 
x < xo ~ 0.1 (or x > xo), and is very steep when x is 
far away from xq. Therefore, the second term in Eq. 
reduces the numerical value of the saddle point when xa 
and Xb are much smaller than xq. As a demonstration, 
let Q = 6 GeV and \/S = 1.8 TeV. Using CTEQ4 par- 
ton distribution and An CD (n / = 5) = 0.202 GeV [0, 
one derives from Eq. (g) that 6$p ~ 1.2 GeV -1 , and 
might conclude that perturbatively resummed prediction 
for Qt distribution at the given values of Q and y/~S is 
not reliable. However, as shown in Fig. [l](a), the inte- 
grand of the 6-integration in Eq. (0) has a nice saddle 
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point at 60 ~ 0.38 GeV -1 , which is within the perturba- 
tive region. This is due to the fact that 0.003 
are very small and the second term in Eq. (||) is negative 
and important. Fig. 0(b) shows that the second term in 
Eq. (a) (dashed line) cancels the first term (solid line) at 
the saddle point bo, which is much smaller than the 6sp 
estimated from Eq. (0). 
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FIG. 1. Integrand of the b- integration in Eq. (Q) at Qt = 
as a function of b for Q — 6 GeV (a) and Q = Mz (c) with 
an arbitrary normalization; and the first (solid) and second 
(dashed) terms in Eq. (Q) as a function of b in (b) and (d) at 
the respective Q. 

In Fig. 0(c) and[j](d), we show the effect of the second 
term in Eq. (^) on the saddle point of Z production at 
the LHC energy. At \/~S = 14 TeV, the dependence 
for Z production improves the value of the saddle point 
from b = 0.24 GeV -1 at VS = 1.8 TeV to 0.13 GeV" 1 , 
in comparison to an estimated 6sp ~ 0.40 GeV" 1 from 
Eq. (J5J) . The narrow width of the 6-distribution shown in 
Fig. |l|(c) also ensures that the 6-integration is dominated 
by b ~ bo- Similarly, we find the same improvements 
on the saddle point and width of the 6-distribution for 
Higgs production at LHC energy [||. In conclusion, the 
6-space resummation formalism for heavy boson produc- 
tion at collider energies should not be very sensitive to 
the nonperturbative input F NP at large b. 

To quantitatively demonstrate the sensitivities on 
the F NP , we reexam the extrapolation, Wcss; de- 
fined in Eq. (|]). We find that using the fitting pa- 
rameters from Rcfs. [Q-[L(| to fix the F NP , the ratio, 
Wcss(b,Q,XA,XB)/W(b,Q,XA,XB), differs from one by 
as much as 20 percents within the perturbative region: 
b < b m ax ~ 0.5 GeV" 1 . That is, the W css {b,Q,x A ,x B ) 
introduces a significant fitting parameter dependence to 
the resummed 6-distribution in the perturbative region. 

In order to separate the perturbative prediction in 
small 6 region from the nonperturbative physics at large 
6, we review the resummation of large logarithms in the 
CSS formalism, and introduce a new functional form for 
the extrapolation. In Rcf. [||, the large logarithms are 
resummed by solving the evolution equation 



2 W(b, Q,xa,x b ) = [K(bfi, a s (n)) 
+ G(Q/^i, a s (n))] W(b, Q, x A , x B ) , 



(7) 



from ln(c 2 /6 2 ) to ln(Q 2 ), and the renormalization group 
equations for K and G from In^c^Jb 2 ) to ln(/x 2 ) and from 
ln(/z 2 ) to ln(Q 2 ), respectively ||. Since the evolution 
equation and the renormalization group equations do not 
include any power corrections, the solution, W in Eq. (||), 
is valid only for 6 < b c with ln(l/6 2 ) ~ 6 2 (or 6 C ~ 
0.75 GeV" 1 ). The choice of b max = 0.5 GeV" 1 in Ref. § 
is consistent with the approximation. 

Taking advantage of our early conclusion that heavy 
boson production at collider energies should not be 
very sensitive to the large 6 region, we extrapolate 
W(b, Q, xa, xb) to the large 6 region without introduc- 
ing the 6» ||. For 6 < b max , the W(b,Q,x A ,x B ) is 
the same as the perturbatively calculated one given in 



Eq. (2 



For 6 > b max , we solve the evolution equation in 



Eq. (^) from ln(c 2 /b max ) to ln(Q 2 ), and solve the renor- 
malization group equations for K and G from ln(c 2 /6 2 ) 
to ln/x 2 and from ln(/x 2 ) to ln(Q 2 ), respectively. Because 
6 > bmax, we add possible power corrections to the renor- 
malization group equations of K and G. We find M 



WQz(b,Q,x A ,x B ) = W(b max ,Q,x A ,x B ) 

X F NP (b,bmax,Q,X A ,X B ) 



for 6 > bmax- Including only the first power correction, 
the F NP has the following functional form M 



\n(F NP ) = - 9l [{b 2 ) a - (b max ) a 



92 



(9) 



where gi, gi and a(< 1) are parameters. Their depen- 
dence on Q, xa, and x B , which are more relevant for 
the low Q 2 Drell-Yan data, are explained in Ref. ||. In 
Eq. (^), the first term corresponds to a direct extrapola- 
tion of the logarithmic contributions to the function K 
to large 6 region. The (6 2 ) a dependence is a result of 
replacing a series of logarithmic dependence on fi 2 in the 
renormalization group equations by {[i 2 ) a ■ The second 
term is a consequence of the first power correction (l//i 2 ) 
to the renormalization group equations. Since the saddle 
point has a small numerical value in 6, high power cor- 
rections to the renormalization group equations of the 
K and G, which are sensitive to the very large 6 re- 
gion, could be neglected. In addition, we neglect a term, 

m (!.r •■ xa, x B )) — \n(W{b 

max 1 W u max j 

X A ,X B )), 

for the hi(F NP ) by assuming that parton distributions 
are saturated when 6 > b ma x- More detailed discussions 
on the functional form of the F NP are given in Ref. [|| . 

We quantitatively test the sensitivities on the F NP by 
studying the dependence on b max , gi, and a. We first 
set gi = (no power corrections) and fix g\ and a in 
Eq. (J)|) by requiring the first and second derivatives of 
the W to be continuous at b = b mnr = 0.5 GeV" 1 . Wc 
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plot our predictions (solid lines) to the Qt distributions 
of Z and W production at Tevatron in Figs. || and [|, 
respectively. In Fig |[ we plot the da/dQx of e + e~ pairs 
as a function of Qt at \/S = 1.8 TeV. The data are 
from CDF Collaboration [O]. Theory curves (Z-only) 
are from Eq. (Q) with the W given in Eqs. (|^) and (j|) 
for b < b max and b > b max , respectively. CTEQ4 par- 
ton distribution and an overall normalization 1.09 are 
used p3[ . In Fig. ||, we plot the da/dQx for W produc- 
tion with the same b max and $2 and without any overall 
normalization ||. The data for W production are from 
DO Collaboration Q. The QCD predictions from the 
6-space resummation formalism are consistent with the 
data for all Qt < Q- Furthermore, we let 172 be a fitting 
parameter for any given value of b max - Although the 
fitting prefers gi ~ 0.4 GeV 2 , the Qt distributions are 
extremely insensitive to the choices of b max and 52 • The 
total x 2 is verv stable for b max e (0.25,0.8) GeV -1 and 
g-2 € (0, 1) GeV 2 . In Figs. || and ||, we also plot the theory 
curves (dashed lines) with g<x = 0.8 GeV 2 (twice of the 
fitting value). Non- vanish 32 gives a small improvement 
to the Qt distributions at small Qt- We also vary the 
value of a in Eq. (^ by requiring only the first derivative 
to be continuous at b = b max , and find equally good the- 
oretical predictions, except very mild oscillations in the 
curves at very large Qt due to the Fourier transform of a 
less smoother 6-distribution. The observed insensitivity 
on b m ax, <?2, and a is a clear evidence that the 6-space 
resummation formalism is not sensitive to the nonpertur- 
bative input at large b for heavy boson production. 



predictions at the LHC energy 
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FIG. 2. Comparison between the 6-space resummed Qt 
distribution and CDF data Jl3[ |. The inset shows the 
Qt < 20 GeV region. 

We notice that the theory curve is below the data at 
large Qt- We believe that it is because we have only the 
leading order contribution to the Y term in Eq. (|l|). At 
large Qt-, the Y term dominates. Similar to the fixed- 
order perturbative calculations, the next-to-leading order 
contribution will enhance the theoretical predictions [|l5| . 
In conclusion, the CSS 6-space resummation formalism 
has a good predictive power for heavy boson production 
at Tevatron energy, and it should provide even better 
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FIG. 3. Comparison between the b-space resummed Qt 
distribution and DO data The inset shows the 

Qt < 20 GeV region. 
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